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Abstract 

This paper analyzes the impact of user mobility in multi-tier heterogeneous networks. We begin 
by obtaining the handoff rate for a mobile user in an irregular cellular network with the access point 
locations modeled as a homogeneous Poisson point process. The received signal-to-interference-ratio 
(SIR) distribution along with a chosen SIR threshold is then used to obtain the probability of coverage. 
To capture potential connection failures due to mobility, we assume that a fraction of handoffs result 
in such failures. Considering a multi-tier network with orthogonal spectrum allocation among tiers and 
the maximum biased average received power as the tier association metric, we derive the probability 
of coverage for two cases: 1) the user is stationary (i.e., handoffs do not occur, or the system is not 
sensitive to handoffs); 2) the user is mobile, and the system is sensitive to handoffs. We derive the 
optimal bias factors to maximize the coverage. We show that when the user is mobile, and the network 
is sensitive to handoffs, both the optimum tier association and the probability of coverage depend on the 
user’s speed; a speed-dependent bias factor can then adjust the tier association to effectively improve 
the coverage, and hence system performance, in a fully-loaded network. 

Index Terms 

Heterogeneous cellular networks, Poisson point processes, range expansion, user association, hand¬ 
off rate, downlink coverage analysis. 
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I. Introduction 

It is now widely accepted that the next generation of wireless communication systems will be 
based on multi-tier heterogeneous wireless access with an all-IP infrastructure Q], 0. With the 
anticipated increase in the number of applications available to a handheld device, e.g., voice, data, 
real-time multimedia, etc. (3), [j4j, mobility management will play an important role in providing 
seamless service to the mobile users moving from one access point (AP) to the other. Client-server 
applications such as email, web browsing, etc., are amenable to short-lived connections and do 
not require sophisticated mobility solutions. Media streams, however, can function normally only 
with a maximum interruption of 50msec; while an interruption of 200msec is still acceptable, 
any longer interruption causes perceptible and unacceptable delays [0. Therefore, on the one 
hand, it is desirable to minimize handoffs between APs to avoid any excessive connection delay 
and call drops; on the other hand, depending on the radio resource allocation scheme, a user 
entering the service area of another AP could cause (and suffer) increased interference levels if 
the handoff is not performed. Therefore, in the move towards multi-tier heterogeneous networks, 
the issues of cell association and handoff must be addressed in an efficient manner to minimize 
the service delay and connection failure for mobile users. 

A. Related Work 

In the context of cellular networks, there is a large body of literature studying the delays caused 
due to handoffs l[5|-0, protocols and effective handoff algorithms JT), |[8ll— ifTOll . and multi-tier 
system design with microcells overlayed by macrocells iflTI . lH2l . If there are enough resources, 
the classic handoff algorithms in a multi-tier network assign users to the lowest tier (e.g., the 
microcells) thereby increasing system capacity [[9]]. To account for mobility, based on an estimated 
sojourn time compared to a threshold, the user is classified as slow or fast, and is assigned to 
the lower or the upper tier respectively l(T3l . [|T4| (these works assume a two-tier network). The 
estimated sojourn time depends on the cell dimensions as well as user information such as the 
point of entry and user trajectory llT5l . Similarly, velocity adaptive algorithms use the mobility 
vector, including both the estimated velocity and the direction, to perform the handoff [[8]. To 
avoid the ping-pong effect due to unnecessary inter-tier handoffs, once the user is classified as 
fast, it remains connected to the upper tier regardless of any changes in its speed lfT6ll . Another 
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alternative is to introduce a dwell-time threshold to take into account the history of the user 
before any handoff decision [|9]|. This technique is based on speed estimation at each cell border. 

Whether the handoff is performed solely by the network controller ifTTIl . or autonomous 
decisions by the user equipment are taken into account IF181 . it is desirable to reduce the 
signalling overhead due to unnecessary or frequent handoffs between the tiers or among the 
APs within one tier. The proposed algorithms mentioned above are mainly applicable in large 
cells. Importantly, the handoff rate, sojourn time or dwell time analysis provided in the literature 
consider deterministic AP locations and a regular grid for the positions of the base stations. 
With the increasing deployment of multi-tier networks, especially small cells in an irregular, 
non-determini Stic manner lfl9ll . handoff analyses for modem heterogeneous networks must now 
take into account the randomness of the AP locations by using random spatial models ll20ll . 

One of the common 2-D spatial models, with the advantage of analytical tractability, is the 
homogeneous Poisson point process (PPP) characterized by the density of nodes, A. In this 
model, the number of nodes in area A is a Poisson random variable with mean A A; the number 
of nodes in disjoint regions are independent random variables, and node locations are mutually 
independent. The accuracy of this model for a two-tier cellular network was examined by Dhillon 
et al. ll2T1 : it was shown that the probability of coverage in a real-world 4G network lies between 
that predicted by the PPP model (pessimistic lower bound) and that by the regular hexagonal 
grid model (optimistic upper bound) with the same AP density. Similar results were reported 
comparing the coverage predictions by the PPP and the square grid model lf22l . Furthermore, the 
PPP model provides a tighter bound at cell edges where the probability of having a dominant 
interferer is closer to that found in an actual 4G network lf22l . Wen et al. Il23ll extended the 
analysis in lITTll . and derived the probability of coverage and the average user throughput in a 
multi-hop multi-tier network. In both works, the maximum received signal-to-interference-plus- 
noise-ratio (SINR) at the user is used as the tier connection metric. The user is then in coverage 
if its received SINR is above the tier’s threshold either through direct connection to an AP 
(infrastructure mode) or through another user (ad hoc mode) |[23l . 

Modeling each tier of a multi-tier network by an independent PPP, Jo et al. Ii24l derived the 
outage probability and the ergodic rate of a multi-tier network with flexible tier association. 
A popular method to achieve this flexibility is by adding, in dB, a bias factor either to the 
average received power ll24l . Il25l or to the instantaneous received SINR 11251 . and use this 
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biased received power or biased SINR as the tier connection metric. If all the APs in the same 
tier have the same transmit power (an assumption common to most works), connecting to the 
AP with the maximum average received power results in the 2-D space tessellation in which 
the AP coverage area is represented by a Voronoi cell. The cell associated with an AP then 
comprises those points of space that are closest to the AP; in turn, the cell size is a continuous 
random variable. An analytical approximation for the Voronoi cell size has been derived in Il26l . 
Applying this formula to a network with the user distribution modelled by an independent PPP, 
the distribution of the AP load in terms of the number of users is derived in a single tier ll27ll . 
and a multi-tier network with flexible tier association ll28l . The authors in ll27ll obtained the 
probability that an AP is inactive and used it to derive the probability of SINR coverage taking 
into account only the received signals from the active APs. Using the proportionally fair model 
where each AP equally divides its available bandwidth amongst its users, the load distribution 
was used to derive the downlink rate distribution at a reference user in a multi-tier network ll28ll . 
and a two-tier network with limited backhaul capacity |[29l . 

The effect of the dominant interferer was studied by Heath el cil. If30l to derive a tractable 
model for the total interference at a specific cell in a heterogeneous network. In this model, 
choosing a fixed-size cell (as the cell under consideration) and a guard radius (hence a guard 
region), the interfering APs form a PPP and consist of all APs lying outside the guard region 
with the nearest one as the dominant interferer. This work allows us to evaluate the performance 
for a “given” cell as apposed to a “typical” cell in the entire network. Clearly, in a multi¬ 
tier network, load distribution through tier association affects the number of users serviced by 
each tier. To deal with the resulting inter-tier and intra-tier interference, two main approaches 
have been proposed: 1) spectrum sharing among tiers but with fractional frequency reuse lf3TI : 2) 
orthogonal spectrum allocation among tiers, thereby eliminating inter-tier interference, e.g., |[32ll . 
If33l . The relation between the allocated fraction of spectrum and the tier association probability 
to maximize the probability that a typical user receives its required rate has been studied in ll33ll . 
It is shown that in the case where an average number of users per AP (only controlled by the 
relative distribution densities and the bias) is assumed, the intuitive solution that the fraction 
of spectrum allocated to each tier should be equal to the tier association probability results in 
essentially zero performance loss. 

In all the works mentioned above in the context of heterogeneous networks modelled by PPP, 


DRAFT 


5 


the performance metric is evaluated at a stationary (but randomly located) user. In this paper, 
we focus on mobility-aware analyses. Since interference management and tier association are 
coupled, handoff analysis becomes important in deriving effective mobility-aware solutions that 
can improve not only the service delay but the overall system capacity through load distribution 
in the network. A comprehensive survey of mobility models and their characteristics can be 
found in Il34l . Il35l and the references therein. They are generally categorized into trace-based 
and synthetic models, and are mainly used to simulate the movement of the mobile users in ad 
hoc networks as realistically as possible. 

To the best of our knowledge, the only work that applied a mobility model in a PPP network, 
and derived an analytical expression for the handoff rate and sojourn time is by Lin et al. Il36l . The 
authors proposed a modified random waypoint (RWP) model] in a single-tier irregular network 
with APs modelled by a PPP, and show that the handover rate is proportional to the square root 
of the AP density. This work defines the handoff rate as the ratio of the average number of cells 
a mobile user traverses to the average transition time (including the pause time). Their analysis 
predicts a slightly higher handover rate and lower sojourn time (overall, a pessimistic prediction) 
compared to an actual 4G network. The handover rate and sojourn time predictions in this work, 
along with the coverage predictions in ll2Tl . |[22ll . imply that the PPP model provides a slightly 
pessimistic but sufficiently accurate analysis while being analytically tractable. 

B. Our Contributions 

Attracted by its applicability and tractability, we use the PPP model for handoff analysis in 
an irregular multi-tier network. Our goal is to analyze the impact of mobility, and to use this 
analysis in deriving effective tier association rules and, hence, load distribution in a multi-tier 
network to minimize this (negative) impact. Similar to |[36l . we consider the handoff rate during 
one movement period. However, we use a different mobility model (as apposed to the modified 
RWP) and a different definition for handoff rate that includes the connection metric. Specifically, 
we define the handoff rate as the probability that the user crosses over to the next cell in one 
movement period. The contributions of our work are as follows: 

'The RWP model (37) is one of the most commonly used mobility models for evaluating the performance of a protocol in ad 
hoc networks. In this model, each node picks a random destination uniformly distributed within an underlying physical space, 
and travels with a speed uniformly chosen from an interval. Upon reaching the destination, the process repeats itself (possibly 
after a random pause time). 
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• First, we derive the handoff rate in a network where the AP locations are modeled by a 
homogeneous PPP. This handoff rate can also be interpreted as the probability that the 
serving AP does not remain the best candidate in one movement period. Based on some 
mild approximations, we simplify this expression; our numerical simulations show that our 
theoretical expression provides reliable results over a broad range of system parameters. 

• Next, we consider coverage in a network with mobile users, assuming that a certain fraction 
of handoffs result in connection failure, i.e., we assume that the outage probability is linearly 
related to the handoff rate derived earlier. We use the signal-to-interference-ratio (SIR) in 
an interference-limited network and a pre-specified SIR threshold to define coverage at a 
reference user. We derive the probability that a mobile user, initially in coverage, remains 
so despite its motion. Interestingly, the results show that the degradation in service - even 
for a fast moving user in a network where most handoffs result in outage - decreases with 
the increase in the SIR threshold (provided that the user is allocated its desired resources). 
This approach provides a tractable model to analyze the impact of mobility; specifically, 
we do not attempt to derive a joint coverage probability distribution across the locations of 
a specific mobile user. 

• Finally, we consider a multi-tier network with orthogonal spectrum allocation among the 
tiers, and derive the coverage probability considering handoffs. We show that the overall net¬ 
work coverage can be improved by adjusting the tier association through the corresponding 
bias factor in a mobility-aware manner, hence, improving system performance in a fully- 
loaded network. This multi-tier coverage expression with handoff can easily be generalized 
to include spectrum sharing across tiers using ll24l . 

The remainder of the paper is organized as follows. In Section HU we describe the PPP model, 
the multi-tier network under consideration, and the tier association metric. The handoff rate for 
a typical user in such a network is derived in Section Hill This is the main result of the paper 
upon which the subsequent sections rely. The probability of coverage in a single tier considering 
a linear cost function for the handoff is derived in Section HV] Section |V] presents the probability 
of coverage in a multi-tier network with and without the handoff analysis. We show how the 
bias factor can optimize the tier association for a mobile user to maximize the probability of 
coverage in a multi-tier network. We conclude the paper in Section fVTl 
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Throughout the paper, we use P(F) to denote the probability of event Y and Ey(-) to denote 
the statistical average with respect to the random variable X. 


II. System Model 

We consider a network comprising K tiers in the downlink. Each tier, indexed by k, is 
characterized as a homogeneous Poisson point process with the tuple {I\. A fc , r k j denoting 
the transmit power, the AP density and the SIR threshold respectively. The tiers are organized 
in increasing order of density i.e., Ai < A 2 ■ ■ ■ < A k- Given the density X k , the number of APs 
belonging to tier k in area A is a Poisson random variable, with mean AX k , which is independent 
of other tiers. Furthermore, all the APs in tier k have the same transmit power P k . 

In our model, transmissions suffer from a constant path loss exponent a. Rayleigh fading with 
unit average power models the channel gain for the received signal from each AP. Since networks 
such as those under consideration here are interference-limited, we ignore thermal noise, and 
use the received SIR as the coverage metric. Finally, for tractability, we ignore shadowing. In 
the coverage analysis, we assume orthogonal spectrum allocation among tiers and a reuse factor 
of one within each tier. Therefore, at a typical user connected to tier k, the set of interfering 
APs include all the APs in tier k except the serving AF0. The results derived in this paper can 
be generalized to allow for spectrum sharing across tiers |[24ll . log-normal shadowing lf38l . and 
any arbitrary fading distribution for the interfering signals ll22l . 

To evaluate the received signal at a typical user, we shift all the point processes so that the user 
under consideration lies at the origin. Despite the shift, each tier still forms a PPP with the same 
density. We use the maximum “biased” average received power as the tier association metric 
where the received power from all the APs of different tiers are multiplied by the corresponding 
bias factor B k , and the user is associated with the tier with the maximum biased average received 
power. Eater in the paper, we will optimize the bias factor to account for mobility. Fet ry denote 
the distance between a typical user and the nearest AP in the jth tier. In this setup, the user 
connects to tier k if: 

k = arg max PjL 0 (rj/r 0 )~ a Bj, ( 1 ) 


2 A reuse factor of greater than one can be accounted for by using a reduced AP density in calculating the interference G2i, 
but will not change the handoff rate. 
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where Bj is the bias factor associated with tier j, L 0 is the path loss at reference distance r 0 , and 
a is the path loss exponent. We use r 0 = lm and L 0 = (An/e)~ 2 where e denotes the wavelength 
at 2GHz. Since all the APs in each tier have the same transmit power and bias factor, the best 
candidate from each tier is the AP closest to the user. Without loss of generality, we set Bi = 1. 
If all the tiers have the same bias factor (or simply, Bj = 1, Vj), the tier association metric 
is the maximum received power. Otherwise, Bj > 1 results in an increased number of users 
connecting to tier j. 

Within the serving tier, the user connects to the nearest AP in that tier. If the user is initially in 
coverage, when it moves, it might fall into the coverage area of another AP at a shorter distance, 
and a handoff occurs. Although the user might be in coverage at both locations, rapid changes 
in the serving AP increases the possibility of connection failure. Fig. Q] shows the scenario 
under consideration, li is the user’s initial location at connection distance r from the serving AP 
denoted by AP S . The user moves a distance v in a unit of time, at angle 9 with respect to the 
direction of the connection, to a new location l 2 at distance R from AP S . This model is mostly 
suitable for a scenario where the user moves at a constant speed or it has small variation such 
that it can be approximated by its mean. 


III. Handoff Rate in a Single Tier Network 


Whether the handoff occurs (Fig. |l(a)| ) or not (Fig. |l(b)[ ) depends on the existence of another 
AP in the circle with the user at the center and radius R. In Fig. |T] C denotes the circle with 
its center at li and radius r; A denotes the circle with its center at l 2 and radius R. The two 
circles intersect in at least one point which is AP S . The excess area swiped by the user moving 
from li to l 2 is denoted by A \ A fl C. For the most part, we assume that the user can move 
in any direction with equal probability. Later in the paper, we show that due to symmetry, the 
probability of handoff for the user moving at angle (27 t — 9) is the same as that for the user 
moving at angle 9 with the direction of the connection. Denoting the corresponding random 
variable as 0, the probability distribution function (PDF) of 0 is then set to be non-zero in 
[0,7r); here, we assume a uniform distribution given by f@(9) = 1/n. 

Let Hk denote the event that a handoff occurs for a user connected to tier k. Throughout the 
paper, we denote the complementary event that a handoff does not occur for the user connected 
to tier k as ///,.. Furthermore, let random variable F4 denote the distance between the user and 
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(a) Scenario where the handoff occurs. 



Fig. 1. Scenario where the user is initially at h, at connection distance r from the serving AP, moving a distance v in the unit 
of time at angle 6 with the direction of the connection; (a) handoff occurs if there is another AP closer than R to the user at 
the new location h', (b) the serving AP remains the closest AP to the user at location h. Hence, handoff does not occur. 


the closest AP in tier k. 

Definition 1. The handoff rate, 'H/fv, A /,.) = P (H k ), denotes the probability of handoff for a 
user connected to an AP belonging to tier k, moving a distance v in a unit of time (speed of v). 

The handoff rate, TLk(v , A^), is a function of v and A& and is given by Theorem [Q 


Theorem 1. Consider a mobile user at connection distance r in a network with access points dis¬ 
tributed according to a homogeneous PPP with density A&. The probability of handoff iP ( If. \ r, 9) 
for the user moving a distance v in a unit of time at angle 9 with respect to the direction of the 
connection is given by: 


P (H k \r,9) 
where R = 



yV 2 + v 2 + 2 rv cos 9. 


7T — 9 + siiT 

Furthermore, 


j / v sin 9 


R 

in such a 


— r 2 (7r 
network, 


— 9) + rv sin 9 J 
the handoff rate 


, ( 2 ) 

for a 
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uniformly distributed 9 is given by: 




1 

1 - - 


PIT POO 

/ / 2n\ k rexp 

1_ 

CN 

^ g 

1 

i _ 

O 

II 

£■» 

O 

II 

-V 

\ 


7r — 9 + sin 


j / n sin 0 


+ r~9 + rn sin 0 


drd6f 


(3) 


Proof: See Appendix |Al ■ 

For the special case 0 = 0 when the user is moving radially away from the serving AP, the 
handoff rate can be written in closed form and is given by the following corollary. 


Corollary 1. The handoff rate for the user moving radially away from the serving AP, i. e., 
fe(9) = 5(9 ) where 5 (■) is the Dirac delta function, is given by: 

9ik(v, A fc ) = 1 - (e~ Xkv2n - 2vTTy/Y k ■ Q(vy/2n\ k )) , (4) 

where Q(x ) = e~ t2 ^ 2 dt. 

Proof: From ([2]). the handoff rate is given by: 


P(H k ) = l-V Rk \P(H k \r,9 = 0) 

no o 

= 1- / e~ Xk < R2 - r2 )f Rk (r)dr 

Jr =C 

/*oo 

= 1 - / 




(0^ _ 


/ r=0 
/•oo 


e -A fc7 r(i?2-r2) . 27rAfcr . e -^A dr 
g-Afc7r(r+w) 2 . 2vr A fc rdr 


(5) 


' r— 0 


g—Afcf 2 7 r _ 2 w 7 r 0 ^ . 


-t- 


I vy/2n\ k 


\[2ff 


'/ 2 dt 


= 1 - (e-^ - 2 n 7 rVAl • Q{vy/^\ k )) , 


where (a) follows from using f? 2 = (r + v) 2 when 9 = 0, and ( b ) follows from the change of 
variable t = \Z2n\ k (r + n), giving the desired result. ■ 

As is clear from ©, in the general case where 9 is uniformly distributed, there is no closed- 
form expression for the handoff rate. However, assuming the user displacement is much smaller 
than the connection distance, v -C R, the handoff rate for the general case can be further 
simplified as derived below. 
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Corollary 2 . The handoff rate 'Hifv. X k ) for a typical mobile user moving a distance v in a unit 
of time in a network where v <C R such that v s |‘‘ 0 ~ 0 and the access points are distributed 
according to a homogeneous PPP with density Xk is given by: 


'Hk(v,X k ) = 1 


71 Je =o L 


1 - 2 b(v, X k , 9)-^ r ne b2 ( v ’ Xk ’ e ' l Q(V2b(v, X k , 6)) exp (- X k v 2 (n - 9 )) d 9, 

( 6 ) 


where b(v, X k , 9) = \ZnX k ^j^-, and a{9) = 2 cos 9{i r — 9) + sin 9. 


Proof: Under the assumption that the movement per unit time is much smaller than the 
connection distance, i.e., vs, jf e ~ 0, the probability of handoff conditioned on r and 9 in © in 
Theorem Q] is simplified to: 


P(H k \r, 9) = 1 — exp — X k R 2 (ti — 9) — r 2 (7r — 9) + rvsin9 


= 1 — exp ( — X k v 2 {y r — 9) + rv (2cos6 ) (7t — 9) + sin (9) 


(7) 


= 1 - exp - X k 


v 2 (tt — 9) + rva{9) 


where a(9) = 2cos9{yi — 9) + sin#. The handoff rate is then given by: 
P(^ fc ) = l-E e [E flh [p(# fc |r,0) 


PTV POO 


= 1 - - I I e- Xk ( v2{n - e)+rva W)f Rk (r)drd9 

71 Jq = o J r= o 

P7T pOO 




77 10=0 Jr =0 


e -A k (vRn-6)+rva{0)) . 2nX k re~ vXkr2 drd9 


=1--/ e 
71 Jo= o 

i r 

77 Jo =0 6 

71 Jo= 0 


- \kV 2 (ir-e) 




2nX k rd:rd9 


I r =0 


e " I t — va{9 ) 



2vr 


dtd# 


-\ k v 2 (-K-e) _ gA/o7r( 2 ^) 2 f i 2 /2 

Jt=V2^^ 

1 - 26(u, A fc , 9)^ r ne b2{ - v ’ Xk ’^Q(V2b(v, X k , 9)) exp (~X k v 2 (7r - 9)) d 9. 


( 8 ) 


Setting b(v, X k , 9) = \/f:X k -f^- and employing the change of variable t = \Z2nX k (r + ppp) 
gives the desired result, and the proof is complete ■ 
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(a) Handoff rate versus user displacement. A k = l/(1000m 2 ). 


(b) Handoff rate versus AP density, v = 5m/s. 


Fig. 2. Handoff rate versus: (a) user displacement in a unit of time v, (b) AP density A k, for both the general case (9 has 
uniform distribution) and the special case of radial movement (9 = 0). 


In Fig. [2l we examine the validity of the expressions for the handoff rate. As seen in the figure, 
since © is an exact expression, the numerical simulations exactly match the analysis for the 
special case where the user is moving radially away from the serving AP. As expected intuitively, 
the handoff rate increases with the increase in the user displacement or the AP density. For the 
general case, the plots are obtained using the approximate expression in ([6]), and the integral 
is obtained numerically. The slight deviation of the analysis from the numerical simulations for 


large v in Fig. |2(a)| or very large AP density in Fig. |2(b)| is the result of user displacement 
becoming comparable to the connection distance which is inconsistent with the approximation. 
However, for reasonable speeds and AP densities, the approximation in © is quite accurate and 
will be used throughout the paper. 

The handoff rate, derived in this section, allows us to incorporate the user mobility in the 
coverage analysis and tier association discussed in the following sections. 
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IV. Coverage Probability with Handoffs 

In If22l . the probability of coverage in an interference-limited, single-tier network with the AP 
locations modelled by a homogeneous PPP, was derived to be: 

1 


P (7 > t) = 


(9) 


1 + p(r, a) J 

where 7 is the received SIR at the user, and r denotes the SIR threshold. For the case where 
both the desired and the interfering signals undergo Rayleigh fading, p(r, a) is given by: 

1 


p( T ,a) = t 


- J2/a 


rdw. 


( 10 ) 


J T -2/ a 1 + U °/ 2 

The expression given in (0) is the probability that a single user is in coverage; it can also be 
interpreted as the fraction of all users in coverage at any given time. Each user connects to the 
AP with the largest signal power. In the analysis of ll22ll . the users are stationary, and there are 
no handoffs. Furthermore, with the focus on SIR, © indicates that the probability of coverage 
is independent of the AP transmit power or density. In other words, increasing the density of 
APs or transmit power will not affect the resulting SIR. Consequently, when handoffs are not 
accounted for, the network capacity increases linearly with the number of APs. However, for 
a mobile user, as shown in Fig. |2(b)[ the handoff rate increases with the AP density which 
negatively affects the system performance. 

The cost associated with handoffs is due to service delays or dropped calls. The higher the 
handoff rate, the higher the chance of degradation in the quality of service. Once we incorporate 
mobility, a fraction of users that initially met the SIR criterion for coverage would experience 
a connection failure due to the handoff. This fraction is determined by the system sensitivity to 
handoffs. Conversely, if the user is in coverage and no handoff occurs, it stays in coverage and 
continues to receive service. To take user mobility into account, we consider a linear function 
that reflects the cost of handoffs. Under this model, the probability of coverage is given by: 


Pc(v, A k , /?, Tfc, a) = P ( 7 fc > Tfc, H k ) + (1 - /?)P ( 7 k > T k , H k ), 


(ID 


where 7*. is the received SIR at the user. The first term is the probability of the joint event that 
the user is in coverage and no handoff occurs. The second term is the probability of the joint 
event that the user is in coverage and handoff occurs penalized by the cost of handoff; here, 
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(3 E [0,1] is the probability of connection failure due to the handoff, i.e., a fraction (3 of handoffs 
result in dropped connections even though the user is in coverage from the SIR point-of-view. 
The coefficient (3, in effect, measures the system sensitivity to handoffs. Its value depends on 
a number of factors, e.g., the radio access technology, the mobility protocol, the protocol’s 
layer of operation and the link speed (EJ—lO- At one extreme, as (3 —> 1, P c (v, A fc , (3, r fc , a) = 
P (o'A: > Tfc, Ilk ), stating that only those users that are initially in coverage and do not undergo a 
handoff maintain their connection, i.e., every handoff results in an outage. At the other extreme, as 
(3 —> 0, the system is not sensitive to handoffs and the expression for the probability of coverage 
reduces to ©, since P (77 > r k ,H k ) + P ( 7 fc > Tk,H k ) = P( 7 fc > T k ). It is worth noting 
that (fill) predicts the coverage probability for a mobile user, already in coverage, immediately 
after knowing the user’s speed and the system sensitivity to handoffs. 

Using the handoff rate derived in Section [TlT] and the overall probability of coverage given 
in CD, we can incorporate the user mobility in the coverage analysis as follows. 

Theorem 2. The probability of coverage P c (v, X k , (3, T k , a) for a typical mobile user moving 
a distance v in a unit of time in a network with access points distributed according to a 
homogeneous PPP with density X k is given by: 


P c (v,\ k ,f3,T k ,a) = 


1 + Pk 


(1 ~( 3 ) + (3 


i r 
71 Je =o 


1 — 2b'\/ne b ' Q(V2b') exp ^ — X k v 2 (n — 9)^jd9 

where b’ = b'(v, X k , 9, r k , a) = and p k = p(r k , a). 

Proof: From (fill) , the probability of coverage conditioned on r and 9 is given by: 


( 12 ) 


P ( 7 k > n, H k \r , 9) + (1 - j3) P (7*. > r fc , H k \r, 9 ), 


(13) 


where 

P ( 7 k > Tfc, H k \r, 9) = P( 7 fc > r k \r) ■ P(H k \r, 9) 

= P(7fc > T-fc|r) ■ (l - P(H k \r,9)) 

^ _ (14) 

= P(7fc > T k \r) - P( 7 fc > r k \r) ■ P(H k \r,9) 

= P(7fc > T k \r) - P( 7 fc > r k ,H k \r,9). 
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Hence, 

P (7 k > r k , H k \r, 9) + (1 - /3)P ( 7 *. > r k , H k \r, 9) 

= P (7 k> r k , H k \r, 9) + (1 - / 3 ) ^P( 7 fc > T fe |r) - P( 7 fc > r k , H k \r , 9)^ (15) 

= (1 -/?)P (7 k > r k \r) +/)P( 7 fc > r k ,H k \r,9). 

The probability of coverage with Rayleigh fading at connection distance r is given by lf22l 
Theorem 2]: 

r^ k >r k \r) = e~^ r2 ^’ a \ (16) 


therefore: 


Rk 


lP(7fc > r k ,H k \r,9) 


(17) 


P c {v,\ k ,P,T k ,a) = E© E /4 P(7fc > T k \r, 0)^ 

= (1 - /3)P(7fc > r k ) + /3E© E 
= i+ P (T k ,a) + E Rk IP(7fc > T k, H k \r, 9) 

In the second term, the probability of the joint event that the user is in coverage and the handoff 
does not occur is given by: 


E 


Rk 

7r /*oo 


IP ( 7 k > Rk: H k ) = E© 

7T Jg =0 J r= 0 

f*7T POO 


(a) - 


7T la— 1 


0=0 Jr =0 


IP (7 k > n, H k \r , 9 ) 

p (7fc > Tk\r) ■ p (H k \r, 9) ■ f Rk (r) drd9 

-tr\ k r 2 p(r k ,a) . & ~^k (v 2 (n-0)+rva(0)) . 271 X k re~ Xk7Tr2 drd9 

e V 


e -\ k v 2 {*-e) e 


71 Je =0 

1/1 


27rAfcrdrd# 


vr V 1 + Pfc 


J r =0 

3 -Afct; 2 (7r-0) 


' 0=0 


1 — 2 b' \Zne b Q(V2b') 


d9, 


(18) 


where p k = p(T k , a) and (a) follows from using the probability of the complementary event 
in CD) given by: 

P [H k \r,9) = 1 — P (H k \r, 9) 

(19) 


= exp - A a 


v 2 (n — 9) + rva(9 ) 


The final two steps are similar to the proof for Corollary [2l employing the change of variable t = 
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va(9) 
~ 2?r 


1+Pfc 


. Finally, using (fl8l) 


y/2n\ k (l + p fe ) (r + 2 ^ 1 +^) ) and setting f/(u, A*, 0, r fc , a) 

in (IT7T) gives the desired result, and the proof is complete. ■ 

The simulation results in Fig. [3] show the effect of mobility on coverage. In both figures, the 
probability of coverage for stationary users is l/(l + p(rfc, a)) = 0.49 for r k = OdB, and a = 3.5. 
For mobile users, on the other hand, the probability of coverage decreases with the increase in 
the handoff rate; as expected, this negative effect is more noticeable when the probability of 
connection failure due to handoffs is significant (i.e., large 3 as in Fig. |3(b)[ ). While the handoff 
rate increases linearly with the user speed in a network with low AP density, it saturates in a 
high density network. 

The plots in Fig. [4] show the probability of coverage versus the SIR threshold r k for a mobile 
user for different AP densities. As expected, although the SIR distribution remains the same 
regardless of the AP density, the probability of coverage in a network with a higher AP density 
decreases due to frequent handoffs. The degradation in coverage not only depends on the network 
sensitivity to handoffs, determined by /?, but also on the SIR threshold. While the probability 
of coverage in a network with a high AP density is lower than that with a low AP density, the 
difference between the two is more noticeable at lower SIR thresholds, or in a network with 
a large f3. It is this interplay between mobility and tier association that leads us to consider 
mobility-aware tier association in the next section. 


V. Mobility-Aware Tier Association 

In the previous section, we derived the probability of coverage for a mobile user in a single¬ 
tier network, and showed how it is affected by the handoff rate. The dependence of the handoff 
rate on the AP density is the rationale for associating fast moving users with the higher tiers 
(with smaller AP densities) to compensate for the potential connection failure due to handoffs. In 
practice, this would be achieved by adjusting the tier’s bias factor. In this section, we first present 
the probability of coverage and the optimum tier association and bias factor for a stationary user 
in a network. Then, as in Section [TV] we incorporate mobility in the coverage analysis assuming 
a linear cost function to reflect the impact of handoffs. Finally, we optimize tier association in 
a multi-tier network to account for the effect of handoffs in maximizing the overall probability 
of coverage for a mobile user . 

In a multi-tier network with the maximum biased average received power as the tier association 
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(a) Probability of coverage versus v; 0 = 0.3. (b) Probability of coverage versus v, 0 = 0.9. 


Fig. 3. Probability of coverage versus user displacement v in a unit of time for different AP densities and Tk = OdB: (a) the 
system is less sensitive to handoffs, 0 = 0.3; (b) the probability of connection failure due to handoffs is large, 0 = 0.9. 



r (dB) 



(a) Probability of coverage versus Tk', 0 = 0.3. 


(b) Probability of coverage versus Tk', 0 = 0.9. 


Fig. 4. Probability of coverage versus SIR threshold Tk for; (a) 0 = 0.3 , (b) 0 = 0.9. v = 15 in both figures. 


metric, the probability that a user connects to tier k is determined by the tier’s AP density, X k , 
transmit power, P k , and bias factor, B k , and has been shown to be ll24l : 


A k = 


A k (P k B k ) 2 / a _ 1 

Ef= 1 Ay (PjBjy/* Ell A APA) 2/a ’ 


( 20 ) 


where P, = jr, B ; j = , X ;] = j 1 . While the transmit power and the AP density are mostly 
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determined by the network infrastructure, adjusting the bias factor can dynamically change the 
user association to different tiers in the network. Let n denote the index of the tier associated to 
the user. Since the user connects to only one tier at a time, the probability that the user connects to 
tier k at connection distance r is given by [24. Lemma 1] P(n = k\r) = riy=i j^k e~ 7TX A p j B A 2/ar2 . 
The probability of coverage for a stationary user in a multi-tier network with the spectrum 
shared across the network is derived in ll24ll . Our focus is on the special case of orthogonal 
spectrum allocation amongst tiers and the corresponding optimum tier association. The optimum 
tier association and the bias factor for the maximum SIR coverage for a single-tier two-RAT^] 
network was derived in ll28l Proposition 1], Generalizing this result to a Ji-tier network with 
orthogonal spectrum allocation across tiers is straight forward and we present it below for use 
later: 


Proposition 1. (a) The probability of coverage for a randomly located user in a multi-tier net¬ 
work with orthogonal spectrum allocation among tiers, and the APs in each tier distributed 
according to a homogeneous PPP with density Xk, is given by: 

K 


Pc S A- 1 


( 21 ) 


^ A k + P( r k, a) 
where Ak is the tier association probability. 

(b) The optimum tier association probability to maximize the SIR coverage is then given by: 

1 /p{Tk,Ot) 


A* - 


( 22 ) 


Efc=l 1 M' r *> Q! ) 

Proof: The proof of Proposition is omitted as it is a straightforward extension of the proof 
in Il28l Proposition 1], ■ 

It is easy to see that when all tiers have the same SIR threshold, i.e., = r V7;\ the maximum 


’RAT refers to radio access technology. Each RAT is allocated a different frequency of operation. 

4 Note that in a multi-tier network with all tiers sharing the same spectrum, the set of interfering APs include all the APs in the 
network except the serving AP in the serving tier. In this setup with the same tier association metric, the overall probability of 
coverage is given by ( 12 1 1 with the term p(jk , a) replaced by X2f=i ^jP^ a Z{rk,a, Bj), where Z(rk, a, Bj ) is given by |24| : 


Z(rk,a, Bj) = T) 


'I, 


(Bjl-rkpn* l + w “ /2 


du 


to account for the interference from the other tiers. For such a case, this expression would replace the corresponding expression 
in the analysis for orthogonal allocation given below. However, it appears that no closed-form expressions are possible and one 
must resort to numerical solutions. 
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probability of coverage is achieved with equal tier association, i.e., A* k 
probability is given by: 

P* = K 
c K + p(r, a) 


1 / K, and the coverage 


(23) 


Importantly, given the tier association probabilities, the optimum bias factors can be found 
uniquely by solving a system of linear equations. 

Define x* = [x 2 , x 3 ,..., x K \ T with x k = B k a , k = 2,3, •••if and x\ — B 1 — 1 for the 
uppermost tier. [ • ] T denotes the transpose operation. From (l20l) . the tier association probability 
for tier k can be rewritten as: 


3 =1 

= 1+ t A \pAbT-BA- 

j=l,j^k 

=* (1 (\Pi L '°) B T = 

j=l,j^k 

~ ~ 2/a 

Setting a,j k = A jPj , the last line in (1241) can be written as: 

K 

(l - A jT 1 ) Xk + ^ cijkXj = -ai k k = 2, 3, • • • if. 

j=2,j^k 


(24) 


(25) 


Given the optimum association probabilities, {A* k } k=l , the optimum vector x* = [x 2 , x 3 ,..., xx] T 
is the unique solution to Ax = b where b = [—ai 2 , — ai 3 ,..., — ai K \ T , and A is given by: 


A = 


(l — Ag ) a 32 

«23 


( 1 -Ar 1 ) 


OR2 

OR3 


0-2 K 0 3 k ••• (1 — A* k *)_ 

Note that this result is not limited to the case of optimum tier association probabilities. Given 
any set of non-zero tier association probabilities, the corresponding bias factors can be found 
solving the system of linear equations above. The following theorem states that the matrix above 
is full-rank and, hence, the relationship between the association probabilities and bias factors is 
one-to-one. 
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Theorem 3. The matrix A is full-rank. 

Proof: See Appendix |B] ■ 

The expressions for the probability of coverage and the tier association provided above do not 
take into account user mobility, the associated handoffs and the connection failure due to such 
handoffs. Using the handoff rate derived in Section HIT] and the linear cost function given in (fTTI) . 
we can generalize the results in Section |TV] to incorporate user mobility in the probability of 
coverage in a multi-tier network as follows: 

Theorem 4 . The probability of coverage P c (v, {Afc}^ =1 , !3. {r k } k=1 , a) for a typical mobile user 
moving a distance v in a unit of time in a multi-tier network with the biased average received 
power as the tier connection metric and the access points of tier k distributed according to a 
homogeneous PPP with density A k is given by: 


Pc(v,{\ k } k= 1 ,P,{r k } k= 1 ,a) = 

K , r, m „ 


A-l 


r[ A k +p(Tk,a) 


(1 - / 3 ) + /3 


71 Je =o L 


1 - 2 b'l 2 Q ( \[2b 


exp 


Afc)n 2 (vr - 6 ») N )d 6 > 


(26) 


where h k = b k (v, X k ,d,r k ,a, A k ) = a) an d o,{6) = 2 cos 6 l ( 7 r — 0) +sin 6 | . 

Proof: Specializing (fl5T) to tier k, we have: 

P ( 7 k >r k ,n = k, H k \r, 6) + (1 - {$)P ( 7 *. > r fc , n = k, H k \r, 6) 

= (1 - /3)P(7 k >r k ,n = k\r) + f3V(y k > r k , n = k , H k \r, 6). 


(27) 


The probability of the joint event that the user connects to tier k, is in coverage and a handoff 
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does not occur is given by: 


P ( 7 k >r k ,n = k, H k ) = E© E Rk P (7*. > T fc , n = k , H k \r, 0) 

-I /*7T /*00 


ft Je =0 7r=o 

^ /*7T /*00 

^ Je=0 Jr=0 


P (7fc > T fc |r) • P(n = fc|r) • P (# fc |u 9) ■ fR k (r) drd0 

e --K\ k r 2 P (T k ,a) _ | e — KXjiPjBj) 2 / 01 ^ j _ e -A fc (w 2 (7r-0)+rm(6i)) _ 27rA/ c re _Afc,rr2 drd0 


0 -\ k v 2 {ix-G) 


ft Je =0 

1 r 


f°° -n\ k r 2 p(r k ,a)+J2f=i>^j{PjBj) 2/o 

/ 6 

7=0 


- • 27rA fc e- A *™ a(0) drd0 


-\ k V 2 (ir- 9 ) 




va(6) )2 

2,r(Afc +Pfe) ^ 27 rAj.rdrd 0 


7T 70=0 
1 ( 1 


^ V Ac 1 + p k y Je= 0 


' r=0 

-A fe i; 2 (7r-0) 


1 - 26' fc 'yTr e b " 2 Q( V26") 


d0, 


(28) 

where p k = p(r k , a). The change of variable t = ^2'ii\ k (A^ 1 + p fc )( r d ■ ... VA —;), and setting 

b k (v, \ k ,0,T k ,a, A k ) = _ 7r i At gives the final expression. Using the sum probability of 

^ V fc “TPfe 

disjoint events, the probability of coverage in a multi-tier network is then given by: 


iC r 


A(u {Afe}f =1 , P, {r k } f =1 , a) = 

fc=i L 


(1 - 0)P(7* > r fc , n = k) + /3P(7 fc >r k ,n = k, H k ) 


K r 


E 


1-/9 A , f> 


L \A k l + Pk J K \A k l + p k ) Jg = 0 


1-1 


0 -\ k v 2 {it-9) 


1 - 2 \J , k ^J?Q(V 2 b 


k) M 


K 


Yl A-l 


A A k + Pk 


(l-P) + P 


1 r - 

ft J9= 0 L 


1 — 2 b”^Ke b * Q(V 2 b'l) exp ( — X k v 2 (n — 0)^0 


(29) 


and the proof is complete. ■ 

As in the single-tier case, if the user is stationary (v = 0), or there is no connection failure 
due to the handoff (J3 = 0), the expression for the overall probability of coverage reduces to the 
expression in Proposition Q] 

To maximize the overall probability of coverage, let the tier association probabilities { A, } f = , 
become the set of optimization variables. The optimization problem with the objective of max- 
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imizing the overall probability of coverage can then be formulated as: 

K 

P* = max [(i - P)fi t i(Ai) + Pfi, 2 (Ai) 

subject to: ^ A t — 1 

2—1 

Aj > 0 i = 1, ■ ■ ■ K, 

where the objective function in (l30l) is the same expression as in (l26l) . Here, /i,i(A) = 1 + ^.p( Tfe a ) 
can easily be shown to be concave with respect to Ai. Due to the complexity of /* : 2 (A;), it is 
not easy to show this to be concave. However, as shown in Fig. 0 numerically 11 ^42 ^ < 0 
over a wide range of system parameters; this suggests that the function is concave for the system 
parameters considered in this paper. Since linear combinations of concave functions (with positive 
coefficients) is concave, we present the following conjecture: 

Conjecture 1. The probability of coverage in a multi-tier network considering handoff derived 
in Theorem\4\is concave with respect to {A}|1 1 - 



Fig. 5. The concavity of the term fi,2(Ai) with respect to Ai for the lower tier, i.e., i = 2 in a two-tier network . 
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The concavity of the objective function, although not leading to a closed-form solution, helps 
us find the optimum tier association probabilities using standard optimization solvers. 

We now show the effect of mobility and more importantly the bias factor on the probability 
of coverage in a two-tier network through numerical simulations. First, for the purpose of 
comparison, the probability of coverage for a stationary user is shown in Fig. [6j {Ai,A 2 } = 
{0.1, l}/(1000m 2 ) and {Pi, P 2 } = {46, 20}dBm denote the tiers’ AP density and transmit power 
respectively. Tier 1 acts as the reference with bias factor B i = 1 and its association probability is 
given by: A x = 1 — A 2 . As expected, in a two-tier network with equal SIR thresholds, the overall 
probability of coverage is maximized when the user connects to each tier with equal probability; 
further, the numerical value of the maximum coverage is independent of the tier AP density 
or transmit power. This, however, is not the case when mobility and handoff cost is taken into 
account. In Fig. [71 we consider a two-tier network specified by {Ax, A 2 } = {0.1,10}/(1000m 2 ), 
{Pi, P 2 } = {46, 20}dBm, f3 = 0.9 and T\ = r 2 = OdB, and obtain the optimum tier association, 
bias factor and the maximum coverage for three different scenarios: 1) “Optimum Bias” is the 
solution to m assuming a concave objective function; 2) “Optimum Bias at v = 0” leads to 
the optimum tier association for a stationary user regardless of its mobility and handoffs derived 
in Proposition [Q 3) “Max-SIR” depicts the scenario where all tiers have the same bias factor, 
Bj = 1 Vj, and the user connects to the tier with the maximum average received power. We also 
compare the obtained results with the optimum solution through a brute force search. As is clear, 
these results suggest that the conjecture stated above is true for the range of network parameters 
considered in this paper. The figure illustrates the importance of accounting for handoffs in 
a multi-tier network. Including the effect of mobility leads to improved coverage by pushing 
fast-moving users to preferentially pick the higher tiers (tiers with lower densities). 

VI. Conclusions 

In this paper, we analyze, and compensate for the impact of user mobility in multi-tier 
heterogeneous networks. We began by deriving the handoff rate for a typical mobile user 
in an irregular multi-tier cellular network. The AP locations in each tier are modelled by a 
homogeneous PPP We showed that the provided analysis matches the numerical simulations 
over a broad range of system parameters, i.e., AP density and user speed. 

The dependence of the handoff rate on AP density, and the associated cost is the main 
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Fig. 6. Probability of coverage in a two-tier network versus A%. Ai = 1 — A 2 , {Ai,A 2 } = {0.1,1}/1000, {Pi,P 2 } = 
{46, 20}dBm and n = T 2 = OdB. The overall probability of coverage is maximized when Ai = A 2 = 0.5. 


motivation in assigning users to different tiers of the network based on their velocity. In a 
multi-tier network where the user connects to an AP of a tier with the maximum biased average 
received power, the probability of tier association has been shown to be a function of AP transmit 
power, tier density and bias factor ||24|. While the transmit power and the AP density are mostly 
determined by the network infrastructure, adjusting the bias factor can change the user association 
to different tiers in the network. We derived the coverage probability with and without accounting 
for mobility for the case of orthogonal spectrum allocation among the tiers. We conjecture that 
optimizing the bias factors is, in fact, a concave problem allowing for efficient solutions. Using 
the optimal bias factors is shown to improve the coverage probability. 


Appendix A 
Proof of Theorem Q] 


Proof: From Fig. |l(b)[ for a typical user initially connected to AP S at distance r and moving 
to the new location l 2 at distance R, a handoff does not occur if there is no other AP closer 
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(a) Overall probability of coverage versus user speed. 



(b) Probability of association to the lower tier versus user speed. 



(c) Bias factor for the lower tier. B 1 = 1. 


Fig. 7. Coverage in a two-tier network with flexible tier association: (a) overall probability of coverage, (b) probabDB^BF 
association to the lower tier and (c) the bias factor for the lower tier. {Ai, A 2 } = {0.1,10}/(1000m 2 ), {Pi, P 2 } = {46, 20}dBm, 
/3 = 0.9 and t\ = T 2 = OdB. 
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v sinG 


Fig. 8. Relationship between r, v and R. 


than R to the user; hence: 


1-P (H k \r,6) 


= p (n(\A\) = 1 1 iv(|^nc|) = l) 
= v(N(\A\Anc\) = o) 


(a) 

= exp 


A fe (|.A\.AnC|)), 


(31) 


where | | denotes the measure of the specified set with \0\— 0, N(-) is the number of APs in the 
specified area, and (a) results from the null probability of a 2-D Poisson process with density 
X k . The handoff rate depends on the amount of the excess area swiped by the user moving from 
1 1 to / 2 given by: 

\A\Anc\=\A\-\Anc\. (32) 


This measure is the same for the user moving at angle 2n — 9 with the direction of the connection. 
Therefore, due to symmetry, we consider 6 being uniformly distributed only in the range of [0, ^ r). 

In plane geometry, the common area between two intersecting circles with radii r and R, 
where the distance between the centers is v, is given by: 


|.Anc| 


From Fig. [8} 



IV ( r + R 


r 2 +v 2 -R 2 \ , p2 __ a —1 ( R 2 +v 2 -r 

2 vr ) ' 11 LUb y 2 vR 


— v ) (r + R + v) {y + r — R) (v 


r + R). 


R 2 = r 2 + v 2 + 2 rv cos 9, 


(33) 


(34) 


and 

r 2 = R 2 + v 2 + 2 Rv cos(7r — 9 + 0) 

= R 2 + v 2 + 2Rv cos(7r — 9 + sin _1 (^^)). 

Using (l34l) and (1351) in (|33T) . the common area is then given by: 


(35) 
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|„4 nC| = r 2 cos- 1 (-cos9) + R 2 cos- 1 (-cosin-O + sin-'i^)) 
- \y/[(r + v) + R] [(r + v) - R][R+ (r - r>)] [R - (r - w)], 


( 36 ) 


where the third term equals: 

1 


V[( r + v) + R] [(r + v) — R][R+(r — v)] [R - (r - u)] 


= - 2 ( [( r + v f ~ R 2 ] [ R 2 - ( r - v ) 2 ] ) 
= — i ^2 rv [1 — cos 9 ] 2n> [1 + cos 9] j 
= —ru sin 0. 


1/2 


1/2 


(37) 


Using the identity cos 1 (— cos (ip)) — n — (p, we obtain: 


\AnC\= r 2 (7r — 9) + i? 2 


0 — sin 


_ 1 ( v sin 9 
~R~ 


— rv sin 9. 


(38) 


Hence, from (|3TI) and (l32l) . the probability of handoff conditioned on r and 9 is given by: 


P(H k \r,9) = 1 -exp(-A fc (|^l\^lnC|)) 


= 1 — exp 


= 1 — exp 


A k I vr R 2 


r 2 {-K — 9) + R 2 [9 — sin 1 ( t,s ™ 6> )) — rv sin 6 


A fc I R 2 


7T — 9 + sin 1 6 ) — r 2 (7r — 9) + rv sin 9 


(39) 


The handoff rate is then written as: 


n k {v,X k ) = P (H k ) = E e E Rk P(H k \r,9) 

1 


r*7T /*oo 

= 1-/ / exp 

71 Je= 0 Jr =0 


-A, R 2 


tv — 9 + sin 


-l 


v sin 9 


R 


—r 2 (ir—9)+rv sin 9 


■f Rk (r) drd# 


(*7T pOO 


= 1-/ / 2nX k r exp 

71 Je= 0 Jr =0 


- A, I R 2 


7i — 9 + sin 


_ 1 f v sin 9 
~R~ 


+ r 2 9 + rv sin 9 


drd9, 

(40) 


where we used the PDF of R k given by f Rk (r) = 27tA k re nXkr2 , and the proof is complete. 
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Proof of Theorem [3] 
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Proof: The proof has two parts. We first show that the determinant of the matrix is given 
by: det A = (—1) A ~ X 1 . The proof is by induction. The statement is true when K = 2, 

1 U=2 A-i 

since the matrix has only one entry, 1 — A^~ l = When K > 2, the coefficient matrix for 

a A-tier network can be written in the form of the block matrix as: 


A 


U V 

w * 


(41) 


where U is a square matrix of size K — 2, V = [a K2 ,a K3 ,... ,a KK ~ i] T is a column vector, 
W = [a 2K , a 3 Ki • • •, O'K-ik] is a row vector and z — 1 — A* K ~ l is a scalar. Using determinant 
of block matrices |[39l , we have: 


det A = (z- l)detU + det(U- VW), 


(42) 


where in the first term det U = (—1) 


K- 2!-Ei 


i =2 

K — l 




by induction. To calculate the second term, 


note that there is a relation between the off-diagonal entries such that a rj = 1 / a Jt and a % k(ikj = 


i f j. Therefore: 


VW 


1 

^23 


«32 

1 


02R-1 «3R-1 


dK-12 
OR-13 


(43) 


Hence, U — VW is the diagonal matrix diag(— A% 1 , —Al \ ■ • •, A* K _ X x ) with det(U — VW) 
= (-i)^- 2 a=2 1 A* 1 . Using m and algebraic manipulation gives the desired result. 

With det A derived above, the numerator 1 — A* > 0, since = 1, and A* G (0,1) 

for 1 = 1,2,... K. Hence, A has a non-zero determinant. Therefore, it is a full-rank matrix with 
rank K — l, and the proof is complete. ■ 
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